Measurement of positronium hyperfine splitting with quantum oscillation 
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Abstract 

' Interference between different energy eigenstates in a quantum system results in an oscillation with a frequency which is propor- 
tional to the difference in energy between the states. Such an oscillation is observable in polarized positronium when it is placed in 
a magnetic field. In order to measure the hyperfine splitting of positronium, we perform the precise measurement of this oscillation 
using a high quality superconducting magnet and fast photon-detectors. A result of 203.324 + 0.039 (stat.) + 0.015( sys.) GHz is 
obtained which is consistent with both theoretical calculations and previous precise measurements. 
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■1. Introduction 

' Positronium (Ps), the bound state of an electron and a 
positron, is the lightest hydrogen-like atom. Since it is a 
purely leptonic system and is thus free from the uncertain- 
ties of hadronic interactions, it is an excellent object for 
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• ^ studying Quantum Electrodynamics (QED), especially for the 
j>^bound state. The two ground states of positronium, the triplet 
(~| State (1^5 1) and the singlet state (I'^o), are known as ortho- 
Q^positronium (o-Ps) and para-positronium (p-Ps), respectively. 
The difference in the energy between o-Ps and p-Ps is called hy- 
perfine splitting (HFS), whose value is about 203 GHz, which is 
I> significantly larger than that for the hydrogen atom (1.4 GHz). 
A theoretical prediction including 0{a^) corrections has been 
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recently obtained using a Non-Relativistic Quantum Electrody- 
namics (NRQED) approach fT]. The result of this calculation 
deviates from the previously measured values iQ^S by a signif- 
icant margin (3.9 cr, 15 ppm). This discrepancy might indicate 
the signal of new physics beyond the Standard Model similar 
to the discrepancy in the measurement of the muon anomalous 
magnetic moment 10] • Both experimental results and the QED 
I— I prediction should be examined again to confirm the discrepancy 



^ ; in a static magnetic neia, tne two states \s - 
■ ■ 1* = 0, m, - 0) mix to the states |h-) and |-) [7]: 
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\+) = c;|i= l,m^ = 0)-HC(> = 0,m, = 0), 
I-) = C^l\s^l,m,,^0) + Cl\s^Q,m-^Q), 
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In the previous experimental results, there may be a systematic error in the 
material effect due to non-thermalized positronium [5J like the o-Ps lifetime 
puzzle (fid. 
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where x 



^hAHFs ' '^HFS is the HFS value, H is the static mag- 
netic field strength, jUb is the Bohr magneton, h is the planck 
constant, and g' = - J4^^) is the g-factor of the elec- 
tron (positron) including the bound state correction |3] . 

The other states 1.? = 1, m, = +1) do not couple with the static 
magnetic field, and thus remain unperturbed. The energy split- 
ting between |-i-) and \s = l,m. = ±1) (the Zeeman splitting) 
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All previous experiments obtained the value of the HFS via 
the formula above by measuring the Zeeman splitting in a mag- 
netic field of a known strength [9]. There are two distinct 
approaches for measuring the Zeeman splitting. The first ap- 
proach, which was proposed in Ref. 1 10], uses an external high 
power light source with a resonant frequency of A^ix (about 
3 GHz in a magnetic field of ^ 0.8 T) to stimulate the transi- 
tion from |s = l,mj - +1) to In this case, the stability 
of the power, frequency of the light source, and the quality of 
the RF cavity are crucial. This approach has been used in many 
previous experiments, for example Mills et al. |2] and Ritter et 
al. 1 3], and has resulted in measurements with accuracies of 
0(1) ppm. 

The second approach proposed by V. G. Baryshevsky et 
al. 101 makes use of the quantum oscillation between \s = 
\,m^ — +1) and Positrons emitted from a fi^ source 

are polarized in the direction of their momentum due to par- 
ity violation in the weak interaction @. Consequently, the re- 



^The polarization ratio P is determined by the initial velocity i' of the 
positron, P = v/c. 
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suiting o-Ps is also highly polarized. This o-Ps is in a super- 
imposed state of 1+) and \s = l,ins = ±1), and thus oscil- 
lates with a frequency Q./2n. This oscillation was observed 
by Baryshevsky ef aZ. in a subsequent experiment 1 1 1]. This 
approach is free from systematic errors due to the high power 
light source and the RF cavity with a high Q- value. Instead, a 
high-performance time-to-digital converter (TDC) is crucial for 
precise measurement of the time spectrum. Those two comple- 
mentary approaches are necessary to understand the discrep- 
ancy in the theoretical prediction and the experimental values 
of the HFS. In 1996, S. Fan et al. (l^ performed an improved 
experiment using this quantum oscillation method, obtaining a 
result of 202.5 + 3.5 GHz. 

In this paper, we greatly improve the accuracy of the mea- 
surement with the quantum oscillation method by using a very 
high quality magnetic field, a fast photon-detection system, and 
high quality TDCs. This method is based on the spin rotation of 
o-Ps (Ps-SR). It is interesting to note that, using the relaxation 
of o-Ps spin, Ps-SR can be used for probing various materials 
in material science research |13]. Since positronium is much 
lighter than a muon, the relaxation processes of positronium 
spin are expected to be much more sensitive to local magnetic 
structures than /i-SrEI. 
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2. Experimental setup 

The top panel of Fig. [T| shows a schematic diagram of the 
experimental setup, while the bottom panel shows a magni- 
fied view of the positronium chamber. A '''^Ge-Ga source, 
whose end point energy is 1.9 MeV, is used. The radioactivity 
is 30 kBq, and is distributed in the active diameter of 9.35 mm. 
A positron passes through a plastic scintillator (NE102) with a 
thickness of 500 //m. The resulting two light pulses are trans- 
mitted in both directions by light guides to two photomultipli- 
ers (PMT; Hamamatsu R5924-70). The positron then stops and 
forms positronium in the silica aerogel target, which is 10 mm 
in diameter and 10 mm long, with a density of 0.11 g/cnr'. The 
surface of the primary grain is made hydrophobic in order to 
avoid the Stark effect from the electric dipole of the hydroxyl 
groups. The plastic scintillator tags the positron emitted along 
the direction of the x-axis, which results in polarized o-Ps along 
the X-axis. The polarization ratio is estimated to be 0.23 by 
Geant4 simulation, in which the geometry, the threshold of the 
plastic scintillator and the velocity distribution of positrons are 
considered. The entire positronium system is contained within 
a chamber evacuated with a rotary pump in order to reduce pick- 
off annihilation 1 14]. 

The magnetic field (z-direction) is provided by a supercon- 
ducting magnet which was originally developed for medical 
NMR use. It has a large bore diameter of 80 cm and an excel- 
lent uniformity of 10 ppm over the volume of the silica aerogel. 
The magnetic field is measured with an NMR magnetometer 
(ECHO-ELECTRONICS, EFM-150HM-AX) which has a caH- 
bration uncertainty of 35 ppm. 



Unfortunately, the lifetime of o-Ps is much shorter than that of a muon. 



Figure 1 : Schematic diagram of the experimental setup. The top panel shows 
the entire experimental setup. The magnetic field direction is along the z-axis, 
while the LaBr3(Ce) scintillators are placed in the yz-plane. The direction of the 
/}* emitted from the ***Ge-Ga source is along the .v-axis. The bold circle is the 
positronium chamber. The bottom panel is a magnified view of the positronium 
chamber, in which the *^Ge-Ga source, the thin plastic scintillator, and the silica 
aerogel are located. The coordinate systems are separately defined in the two 
panels. 



The produced o-Ps decays into three gamma rays, which 
are then detected by six LaBr3 crystals of 1.5 inches in di- 
ameter and 2 inches long. PMTs (Hamamatsu R5924-70) are 
attached to these crystals. The LaBrs detectors are located 
at (0,0) = (|,-|),(^,-|),(^,|),(|,|),(|,-|),(|,|), 

where 6 - arccos( .i ) and d) = arctan(-). The detectors 

are referred to as LaBrj - 1 to LaBrs - 6 (see Fig.[TJ. The quan- 
tum oscillation modulates the angular distribution of the three 
gamma rays emitted from the o-Ps decay, which leads to beats 
in the decay curve of o-Ps. Unlike the muon precession, in 
which the emission direction of jt/ — > e rotates, this oscillation 
changes its angular distribution as a vibration in the yz-plane. 
This is a unique property of a spin-1 system. The LaBr3 - 1 and 
LaBr3 - 3 detector pair observes the oscillation with the same 
phase, while the LaBr3 - 2 and LaBr3 - 4 detector pair observes 
the inverse phase. The LaBr3 - 5 and LaBr3 - 6 detectors ob- 
serve the exponential decay curve without the oscillation. 

Fig. |2] shows the energy spectrum measured with one of 
the LaBr3 crystals. We note that a good energy resolution of 
4.0 % (FWHM) at 511 keV is obtained. The time resolutions 
of the LaBr3 detectors for the 511 ke V gamma peak and the 
positron tagging plastic scintillator are 200 ps (FWHM) and 
3.8 ns (FWHM), respectively. These values are obtained with 
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Figure 2: Energy spectrum of a LaBrj crystal measured in tlie magnetic field of 
100 mT. 

the PMTs located in a magnetic field of 100 mT. 

Data acquisition starts with the rate of 740 Hz when the plas- 
tic scintillator signal is coincident within -50 ns to 1650 ns with 
at least one of the LaBrs signals. The time f = is defined as the 
timing of the plastic scintillator pulse. A charge ADC (CAEN 
C1205) is used to measure the energy information of the LaBrs 
crystals, while another charge ADC (REPIC RPC-022) is used 
to measure both the baseline information of the LaBr3 crys- 
tals and the energy information of the plastic scintillator. The 
baseline information is measured just before the gamma ray ar- 
rives at the LaBr3 crystals in order to remove pile-up events. 
The time differences between the plastic scintillator and LaBr3 
scintillators are measured by direct clock TDCs (5 GHz: time 
resolution of 200 ps). These TDCs have excellent integral and 
differential linearities. 

Separate measurements have been made for 5 different mag- 
netic field strengths: mT, 100 mT, 118 mT, 135 mT and 
138 mT. Also, both +x (up) and -x (down) polarization mea- 
surements were performed by turning the positronium chamber 
upside down. The expected time periods of the oscillation are 
about 26 ns and 14 ns for 100 mT and 138 mT, respectively. The 
period of each run was about 3 days. In total, 1.4 x 10*^ events 
were recorded in about 22 days of data acquisition period. The 
energy and timing spectra are calibrated every one hour using 
the 511 keV and pedestal peaks. 

3. Analysis 

The following event selections are applied in order to obtain 
a clean time spectrum: 

1. In order to remove pile-up events, the fluctuation of the 
baseline of the LaBr3 is required to be smaller than 
3cr (where cr is the noise level). 

2. The events in which more than two LaBr3 crystals are 
hit simultaneously are disregarded. This helps to reduce 
the accidental contribution since accidental events have a 
back-to-back topology. 

3. In order to obtain a good time resolution, the energy de- 
posited in the LaBr3 is required to be larger than 100 keV. 



Table 1 : An example of the fitted result using the separate fitting method for 
100 mT (up). 



Parameter 


Fitted Value 


Ai 


0.095 + 0.011 


By 


0.078 + 0.012 


Ci 


0.0096 + 0.0003 




0.00733 + 0.00001 


Oi (rad) 


0.18 + 0.03 


n (ns"') 


0.00733 + 0.00012 


72 (ns"') 


0.00975 + 0.00024 


Q. (rad/ns) 


0.24573 ± 0.00017 




1.00 (ndf=8370) 



Since the multiple hit events are removed, six statistically in- 
dependent time spectra are obtained. We then fit the spectra 
using two distinct methods: the separate fitting method and the 
subtracting method. 

3.1. Separate fitting method 

Fig. [3] shows the timing spectra with the best fitted result us- 
ing the separate fitting method. In this method, the six time 
spectra are simultaneously fitted in the range 50 ns to 1450 ns 
with the following functions: 

/„(f) = A„e-y^' + B„e-y^-' 

y\+y2 . 

+ C„e 2 X sin(Qr + 6„) 

+ D„ (6) 

forn = 1,2,3,4 
/„(f) = A„e-y" + B„e-^'-' 

+ D„ (7) 

forn = 5,6 

where A and B are proportionality constants for the decay 
curves, C is the oscillation amplitude, and D is a constant for 
accidental hits (n denotes the index of the LaBr3 detector). The 
two decay rates yi and 72, and the angular frequency of the 
oscillation D. are common variables among the six functions, 
while the others are kept free, yi and 72 are decay rates for 
\s = l,m. = +1) and |+), respectively. These rates include the 
effect of pick-off annihilation. Q. is the Zeeman splitting fre- 
quency defined in Eq. (|5]l. 

The results of the fits for LaBr3 - 1 are listed in Table[T]for the 
100 mT case. All fitted variables converged as shown in the ta- 
ble, in which a reasonable ;if^/ndf of 1.00 is obtained. The fitted 
lifetime values I/71, I/72 are 136.4+2.2 ns and 102.5+2.5 ns, 
which are consistent with the lifetime (1.? = l,m- - +1)) mea- 
sured in aerogel lITill and the calculated value (|+)) in a mag- 
netic field of 100 mT, respectively. A fitted time period In/Q. 
of 25.57 + 0.02 ns (700 ppm) is obtained for this run. 

The same procedure is applied for the runs with different 
magnetic field strengths and different polarizations. The result- 
ing ;^'^/ndf values are always found to be less than 1.03. Ahfs 
can be calculated from the fitted Q.The obtained Ahfs are listed 
in Table |2]for the various magnetic field strengths. 
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Figure 3: The timing spectra obtained in the separate fitting method for 
100 mT (top) and 135 mT (bottom). In both panels, the data points are plotted 
with error bars, while the solid lines show the best fitted results. The opposite 
phase spectra are superimposed in both panels. The polarization direction of /J^ 
is upwards. 



The 118 mT (up) measurement is also performed using a dif- 
ferent TDC clock of 8 GHz, in which consistent results are ob- 
tained. This is an important check for the TDC, since the time 
spectrum is crucial for this experiment. 

The Ahfs values obtained at the various magnetic field 
strengths are consistent with each other. The combined value 
is 203.336 ± 0.048 (stat.) GHz. 

3.2. Subtracting method 

In this method, the sum of the LaBr3 - 2 and LaBr3 - 4 spec- 
tra are subtracted from the sum of the LaBra - 1 and LaBrs - 3 
spectra. Ideally, this would cancel the exponential components 
in the spectra, leaving only the oscillating component. How- 
ever, the acceptances of the LaBr3 detectors are not exactly the 
same, which results in small exponential components remain- 
ing after the subtraction. Although these components are thus 
still included in the fits, the upshot is that the oscillating com- 
ponent is greatly enhanced. Furthermore, the cancellation of 
the prompt peak means that the fitting region can be extended 
closer to zero. Therefore, this method has a smaller statistical 
error in the fits. Fig. |4] shows examples of the subtracted time 




Figure 4: Time spectra obtained in the subtracting method for 100 mT (top) 
and 135 mT (bottom). In both panels, the data points are plotted with error 
bars, while the solid lines show the best fitted results. 



spectra with the best fits superimposed. The fitting region is set 
from 16 ns to 1416 ns. The fitting function is given by 



D, 



' X sin(Q/ + 0) 



(8) 



where yi, 72, and Q are the same as the ones in the separate 
fitting method. The two exponential components with propor- 
tionality constants A and B are for the remnant decay curves, 
while the component with constant C is the oscillation contri- 
bution. The amplitudes A, B and D are expected to be small. 

The fitted results are listed in Table |3] for the 100 mT case. 
The coeflicients A and B are consistent with zero and D is also 
much smaller than C, which means that the cancellation works 
well. A fitted time period In/D. of 25.59 ± 0.01 ns (590 ppm) is 
obtained for this run. 

The obtained Ahfs are listed in Table |4] for the various mag- 
netic field strengths. The Ahfs values obtained at the various 
magnetic field strengths are consistent with each other The 
combined value is 203.324 + 0.039 (stat.) GHz. 
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Table 2: Summary of the HFS values obtained using the separate fitting method, (up) and (down) denote the direction of the y3+ with respect to the .v-axis. (TDC) 
denotes the run in which a different TDC clock was used. Runs in the magnetic field of 1 1 8 mT have larger errors because the runs are performed for shorter periods. 



Run 


Magnetic Field (mT) 


Q. (rad/ns) 


HFS (GHz) 


Events 


lOOmT (down) 


100.592 


0.24588 ± 0.00015 


203.04 + 0.12 


2.7 X 


10^ 


lOOmT (up) 


100.594 


0.24573 + 0.00017 


203.17 + 0.14 


2.1 X 


10« 


1 1 8mT (down) 


118.824 


0.34242 + 0.00032 


203.42 + 0.19 


4.9 X 


10^ 


118mT(up) 


118.826 


0.34289 + 0.00036 


203.14 + 0.21 


6.5 X 


10^ 


118mT(up,TDC) 


118.826 


0.34207 + 0.00039 


203.63 + 0.23 


6.5 X 


10^ 


135mT (down) 


134.805 


0.44034 + 0.00025 


203.58 + 0.11 


1.9 X 


10« 


135mT (up) 


134.807 


0.44104 + 0.00030 


203.26 + 0.14 


2.1 X 


10» 


138mT (down) 


138.326 


0.46394 + 0.00031 


203.45 + 0.14 


9.7 X 


10^ 


138mT(up) 


138.330 


0.46414 + 0.00027 


203.37 + 0.12 


2.0 X 


10« 



Table 4: Summary of the HFS values obtained using the subtracting method, (up) and (down) denote the direction of the with respect to the .v-axis. (TDC) 
denotes the run in which a different TDC clock was used. Runs in the magnetic field of 118 mT have larger errors because the runs are performed for shorter periods. 



Run 


Magnetic Field (mT) 


O (rad/ns) 


HFS (GHz) 


Events 


lOOmT (down) 


100.592 


0.24579 + 0.00012 


203.11 +0.10 


2.7 X 


10*^ 


lOOmT (up) 


100.594 


0.24558 + 0.00014 


203.30 + 0.12 


2.1 X 


10« 


1 18mT (down) 


118.824 


0.34248 + 0.00027 


203.38 + 0.16 


4.9 X 


10^ 


118mT(up) 


118.826 


0.34289 + 0.00028 


203.15 + 0.17 


6.5 X 


10^ 


118mT(up,TDC) 


118.826 


0.34221+0.00031 


203.55 + 0.18 


6.5 X 


10^ 


135mT (down) 


134.805 


0.44061+0.00020 


203.46 + 0.09 


1.9 X 


10« 


135mT (up) 


134.807 


0.44103 + 0.00023 


203.27 + 0.11 


2.1 X 


10^ 


138mT (down) 


138.326 


0.46401 + 0.00024 


203.42 + 0.11 


9.7 X 


10^ 


138mT(up) 


138.330 


0.46422 + 0.00022 


203.34 + 0.09 


2.0 X 


10^ 



over the range of the aerogel. This effect is estimated at 
10 ppm. 

3. The accuracy of the TDC is determined by that of the 
clock (Hittite HMC-T2000), which is better than 2 ppm. 
The effects of differential and integral non-linearities in 
the TDC are negligible. 

4. The fitting region dependence is negligible as long as the 
fitting start time is later than 50 ns for the separate fitting 
method, and 16 ns for the subtracting method. 

5. In the subtracting method, remnant exponential compo- 
nents could affect the HFS value in fitting procedures. This 
amount is estimated by altering the scale factor for each 
histogram by 5%, which result in the shift of 10 ppm in 
the HFS value. In the subtracting method, weighted aver- 
age of the remnant exponential components over all runs 
is only 0.7%, which causes a negligible effect on the HFS 
value. 

6. An electric field produced by aerogel grains varies the 
overlap of the wave function of positronium, which results 
in a shift in the HFS (Stark effect). Since a hydropho- 
bic silica aerogel (hydroxyl groups are replaced with tri- 
methyl-silyl groups) is used in this experiment, the effect 
of Stark shift is rather small. The amount of Stark effect is 
estimated as follow |15]: 

The area density of silanol groups of the aerogel (cr = 
0.44 nm"^) is measured in Ref. fl5j. Although the den- 
sity of aerogel is different, the primary grain in the both 
aerogels is the same. Furthermore, the same chemical pro- 
cesses to replace silanol group is applied. The average 



Table 3: An example of the fitted result using the .subtracting method for 
100 mT (up). 



Parameter 


Fitted Value 


A 


0.0056 + 0.0042 


B 


-0.0015 + 0.0043 


C 


0.0377 + 0.0008 


D 


0.00013 + 0.00002 


71 (ns"') 


0.0077 + 0.0010 


72 (ns"') 


0.0094 + 0.0010 


Q (rad/ns) 


0.24558 + 0.00014 


(rad) 


0.16 + 0.02 


;r'/ncif 


1.00 (ndf=1392) 



4. Discussion and result 

4.1. Systematic errors 

The systematic errors are summarized below: 

1. Varying the frequency sweep range of the NMR magne- 
tometer results in slightly different readings. The uncer- 
tainty in the magnetometer calibration is estimated from 
this deviation (35 ppm). 

2. The non-uniformity of the magnetic field results in the fol- 
lowing two effects: (1) The oscillations in the time spectra 
become smeared. This effect is already taken into account 
in the fitted results listed in Table 2 and 4. (2) There may 
be a difference between the value of the magnetic field 
strength as measured by the NMR and the actual values 
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electric field which affects positronium during its flight can 
be calculated. 



|£|2 - 



1 \2 2 

1 \ mrp 



= l.Ox 10'^ V^/m^ 



(9) 



(10) 

where p — 1.7 x 10""* esu-cm is electric dipole moment 
of hydroxyl groups [ 16], L = 130 nm is the mean distance 
between the grains, and ops - 0.106 nm is the Bohr ra- 
dius of positronium. The amount of the shift in the HFS is 
then 



01 p i 
IB 



^HFS 
HFS 



-248 



-^0 



(11) 

(12) 



= -lOppm, 

where £0 = m^e^ jfi^ - 5.14 x 10** V/cm. This model 
reproduces the experimental values in Ref. |18] within a 
factor 2. This Stark effect is treated as the systematic error 
in this analysis. 

Since Ahfs depends on the magnetic field squared, the sys- 
tematic errors due to the magnetic field uncertainty are doubled. 
The systematic errors are combined by summing in quadrature. 

4.2. Discussion 
Ahfs - 

203.336 ± 0.048 (stat.) ± 0.015 (sys.) GHz {separate fitting) 
203.324 ± 0.039 (stat.) ± 0.015 (sys.) GHz {subtracting) 

We note that the results of the two different fitting methods 
are consistent. We use the more accurate subtracting method 
value as the final result. The accuracy is 200 ppm, which is 
an improvement by a factor 90 over the previous experiment 
which used the oscillation method il2|. This result is consistent 
with both the theoretical calculation flj and the previous more 
precise experimental values which directly measure the Zeeman 
transition 10, H. 

In order to observe the relaxation of positronium spin (Ps- 
SR), the oscillation amplitude is fitted as a function of time. 
The result is, however, consistent with a constant. A better ac- 
curacy and a higher density target are necessary to observe the 
relaxation. 

The accuracy of the measured HFS value in this study is 
200 ppm. The following five points can be improved in order to 
achieve a better accuracy: 

1 . Increase the total run time by a factor 20 (about 1 .5 years). 

2. Increase the intensity of the radioactive source by a factor 
3. The dead time of the DAQ would still be acceptable. 

3. Increase the coverage of the photon detectors by a factor 
3. Fine segmentation is still necessary. 

4. The absolute calibration of the NMR and the uniformity 
of the magnetic field can both easily be improved to 
0(1) ppm. 

5. The Stark effect can be estimated by changing the density 
of the aerogel. 



The result of these improvements would be an increase in 
statistics by a factor 180 and an improvement of the final accu- 
racy to about 15 ppm. 
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